Critical Number of Flavours in QED 
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We demonstrate that in unquenched quantum electrodynamics (QED), chiral symmetry breaking 
ceases to exist above a critical number of fermion flavours Nj. This is a necessary and sufficient 
consequence of the fact that there exists a critical value of electromagnetic coupling a beyond which 
dynamical mass generation gets triggered. We employ a multiplicatively renormalizable photon 
propagator involving leading logarithms to all orders in ot to illustrate this. We study the flavour 
and coupling dependence of the dynamically generated mass analytically as well as numerically. We 
also derive the scaling laws for the dynamical mass as a function of a and Nf. Up to a multiplicative 
constant, these scaling laws are related through (a,Qc) {1/ N f ,1/ N'j). Calculation of the mass 
anomalous dimension 7^ shows that it is always greater than its value in the quenched case. We 
also evaluate the /3-function. The criticality plane is drawn in the (a, Nf) phase space which clearly 
depicts how larger Nj is required to restore chiral symmetry for an increasing interaction strength. 
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The /3-function determines the running of the coupling 
constant. In quantum chromodynamics (QCD), its evo- 
lution, both in the ultraviolet and infrared, is crucially 
influenced by the number of light quark flavours. Virtual 
quarks and gluons contribute to its perturbative tail in a 
diametrically opposed manner and the value of N f deter- 
mines which would be the dominant effect. QGD exhibits 
asymptotic freedom because Nf happens to be less than 
a critical value of N'j^ ~ 16.5. Lattice studies in the 
infrared indicate that just below this value, chiral sym- 
metry remains unbroken and colour degrees of freedom 
are unconfined lli. Below this conformal window, for an 



function renormalization in the bare vertex approxima- 
tion was obtained in In this article, we demonstrate 
that the unquenched QED also has a critical number 
of flavours N'j above which chiral symmetry is restored. 
The starting point is the SDE for the electron propagator 

S-Hp) = S^\p) + ie^ [ d^krS{k)T%k,p)/\^M , (1) 



where q = k — p, e is the electromagnetic coupling and 
Sq^{p) = is the inverse bare propagator for mass- 
less electrons. We parameterize the full propagator 
S{p) in terms of the electron wave function renormal- 



8 < A^^^ < ~ 12 , the evolution of the beta function in the ization and the mass function M(p2) as S{p) 



infrared is such that QCD enters the phase of dynamical 
mass generation (DMG) as well as confinement. QCD is 
not the only gauge theory where infrared dynamics re- 
sponds to the number of fermion flavours in such a dra- 
matic fashion. It has been established that QED3 possess 
a critical number of flavours N^^ associated with the si- 
multaneous emergence of dynamical masses and confine- 
ment if the electron wave function renormalization, pho- 
ton vacuum polarization and electron photon vertex are 
homogeneous functions at the infrared momenta, 0, Q . 
See [3 and [3| for some original works on DMG in QCD 
and QED3, respectively, through the Schwinger-Dyson 
equations (SDEs). 

One ponders if such criticality also characterizes other 
gauge theories in a similar manner. In this paper, 
we study the flavour dependence of DMG in QED4 or 
QED, For large a, it is known to exhibit chiral 

symmetry breaking in the one loop approximation of the 
photon propagator, 0|. A consistent solution for coupled 
equations for the fermion mass function and photon wave 



— M(p^)). A^y(g) is the full photon propaga- 
tor which can be conveniently written as 



A^^(g) = 5— 1 .g^^ 



g2 



(2) 



where ^ is the covariant gauge parameter such that ^ = 
corresponds to the Landau gauge. G{q^) is the photon 
renormalization function. The full electron photon ver- 
tex is represented by V^{k^p). The form of the full ver- 
tex is tightly constrained by various key properties of the 
gauge theory, e.g., multiplicative renormalizability of 
the fermion and the gauge boson propagators, [lol □J, 
perturbation theory, [l3|, t he requirements of gauge in- 
variance/covariance, |13l4l7| and of course, observed phe- 
nomenology, [3 ■ The most general decomposition of this 
vertex in terms of its longitudinal and transverse compo- 
nents is 



r^(fc,p) = ^A,(fc,p)Lf(fc,p) + ^r,(fc,p)7^f(fc,p) ,(3) 
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TABLE I; We tabulate analytical (indicated with superscript 
A) and numerical values (indicated with superscript A'') of Qc 
for different values of Nf, 



and Nf for different values of a. 



where 
and = 



7^ 



til/ 



{k 



^2 



p)^, where a'^" = 



[7'',7'']/2. The co- 



efficients Xi are determined through the Ward-Takahashi 
identity relating the electron propagator with tire elec- 
tron photon vertex |2^ : 
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(4) 



and X4{k,p) = 0. A simple choice of the transverse coeffi- 
cients which, combined with the longitudinal component, 
renders both F{p^) the G{q^) multiplicatively renormal- 
izable for massles fermions, has been constructed only 
recently [ll|. The longitudinal coefficient Ai plays a cru- 
cial role in ensuring the correct leading logarithms are 
summed up for the photon wave function renormaliza- 
tion. Similarly A2 dictates the multiplicative renormaliz- 
ability of the photon pro pag ator. Using this information, 
the ansatz proposed in [il| makes use of the following 
four transverse basis vectors as suggested by Ball and 
Chin, Hi : 

T^{k,p) = p^k-q-k^p-q, 



T^{k,p) = -rik'-p^) + ik+prfi, 

T^{k,p) = -^^k'-p'^axu + k^ji-p^'l^ 



(5) 



The corresponding coefficients are chosen to depend 
upon F{p'^) in the following simple manner : 
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The coefficient A3, which enters into the description of 
massive fermions, is irrelevant to the power law behaviour 
of both F(p2) and G{q^). However, it is intimately re- 
lated to the value of the anomalous dimension 7^ for 
the fermion mass function. In the quenched theory, the 
ultraviolet behaviour of M{p^) can be expressed as 

M(p2) ^ (p2)7™/2-l (7) 

in the deep Euclidean region. At criticality, the mass 
function behaves as Eq. (O at all momenta. If the trans- 
verse vertex vanishes in the Landau gauge, 7,„ = 1.058, 
see e.g. [l^. However, Holdom and Mahanta [2^, using 
the arguments based on the Cornwall- Jackiw-Tomboulis 
(CJT) effective potential technique, have shown that 7^ 
is strictly equal to 1. The importance and usefulness of 
employing the bare vertex was also stressed in (2l| . If it 
were true that "f„i ~ 1, this would suggest that there is 
a necessary piece in the transverse part of the effective 
vertex which does not vanish in the Landau gauge. Com- 
plete calculation of the fermion-boson vertex at the one 
loop in arbitrary gauge and dimensions, reveals that 
the transverse part of the vertex indeed does not vanish 
in tire Landau gauge in any space-time dimensions. This 
fact may possibly favour Holdom's arguments. Thus it 
may well be that the non-zero transverse piece in the 
Landau gauge cancels out the A3 piece of the longitudinal 
component in the equation for the mass function. Con- 
sidering this argument, one such vertex was constructed 
in Following suit, consider tire following full vertex 

r^(fc,p) = r''scik,p) + r^Kpik,p) + r'X{k,p) . (8) 

As the subscripts indicate, T'^u{k,p) is the longitudinal 
Ball-Chiu vertex, defined by Eq. (U), T^p{k,p) is the pro- 
posal by Kizilersu and Pennington, Eq. ([5]), and r^(fc,p) 
is the additional transverse piece which minimally en- 
sures 7m = 1 in the quenched case. With this choice of 
the full vertex, we obtain, in the massless limit 



Pip') = ( ^ 



P 



(9) 



where v — a^/(47r), s — aNf/{3TT), a — e2/(47r) and 
Nf is the number of massless fermion flavours. All main 
conclusions are robust under different truncations, e.g., 
for 1-loop logarithmic photon propagator and for a re- 
summation of the propagators beyond leading logs. Near 
criticality, where the generated masses are small, one can 
assume that the power law solutions for the propaga- 
tors capture the correct description of chiral symmetry 
breaking. We choose to study the resulting equation for 
the mass function in the convenient Landau gauge. Re- 
sults for any other gauge can be derived by applying the 
Landau-Khalatnikov-Fradkin transformations [5,' 17, 24|. 
The usual simplifying assumption G(q^) ~ G{k'^) for 
k"^ > p2 and G{q'^) = G{p'^) for p^ > fc2 allows the ana- 
lytical treatment of tire linearized equation for the mass 
function : 
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FIG. 1: The mass function for different values of a. 



where g{q^) ~ 3aG{q'^)/ (Att) and we have introduced 
ultraviolet cut-ofT A^. The infrared cut-off mimics 
the M'^{k'^) term in the denominator which has been 
dropped off. It is already known that for the one loop 
photon propagator, there exists a critical coupling etc 
above which masses are dynamically generated. One 
can formally demonstrate that the existence of critical 
coupling implies the existence of a critical number of 
flavours above which chiral symmetry is restored. Note 
that in Eq. G{q^) = G{q^,aNf ). Instead of work- 

ing with the variables {a,Nf), if we define a' = aNf, 
we could equally work with {a',Nf). In such case, 
g{q^) = 3a'G(g2)/(4Ar/7r) with G(g2) = G{q^,a'). If we 
hold a' constant, the effective coupling is l/Nf. There- 
fore, the presence of an ac implies the existence of an 
N^. The critical behaviour should thus translate as 
{a,ac) {l/Nf,l/Nj). Moreover, if there is no crit- 
ical a, there will be no critical Nf. We shall demon- 
strate this explicitly for the multiplicatively renormaliz- 
able photon propagator. Let us make the change of vari- 
ables X = A^/p^ and convert the integral equation ([TU]) 
into a second order differential equation 



xHl"{x) + sxM'ix) + — (1 - s)^^ = 



fll) 



with the following infrared and ultraviolet boundary con- 
ditions respectively 

M(l) =A/'(l)/(l-s) , M'{A^/m^) = Q. (12) 



Eq. pT|) can be converted into a Bessel equation through 
Lommel transformations : z ~ Bx'^ , W = x~^ M. Thus 
we work with W{z) instead of M{x). The corresponding 
equation is 



FIG. 2: The mass function for a fixed a and varying Nf. 
As Nf is reduced, the mass function drops significantly for 
increasingly small variations in Nf, suggesting the existence 
of a critical number of flavours. 



where 7 = -s/2, A = (1 - s)/2, A = (1 - s)/s and B = 
\/'ia{l — s)/{tts'^). Moreover, we have assumed s < 1. 
Eq. has the solution 



W{z) = GiJa{z) + C2Ya{z) 



(14) 



where Ja{z) and Ya{z) are the Bessel functions of the 
first and the second kind respectively. The boundary 
conditions get translated as 



W'iz) = '^—^ W{z) 



7B 



■^2 \ 7 

aW{z)+jB[^] W'{z) = 



.(15) 



These conditions allow us to find the constants Ci and 
C2, and the equation for the mass m : 

2XJa{z) + iB (AVm2)^ [Ja-i{z) - Ja+i{z)] 



z^W"{z) + zW'{z) + {z^ ~ A^)W{z) = 



(13) 



2\Ya{z)+^B {Aym^y [Ya-i{^) - Ya+i{z)] ^^^^^A^y 

_ jB [JA-ijB) - JA+ijB)] - 2(1 - s - X)Ja{B) 

7S [Ya-i{B) - Ya+i{B)] - 2(1 - s - X)Ya{B) ' ^ ^ 

Critical a or Nf can be obtained from this equation by 
requiring it to hold true for A — > 00. This implies finding 
the zeros of the equation 

7B [Ja-i{B) - Ja+i{B)] - 2(1 - s - \)Ja{B) = 0.(17) 

For various values of Nf , analytical values of ac have been 
tabulated in the second column of Table. HI Similarly, for 
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The scaling with a for N,=2,0 



The scaling with N, for a = 2.5 
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FIG. 3: The scaling law for the coupling a. Numerical 
solution (solid line) is compared with analytical prediction 
(dashed line) of the linearized equation. 
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TABLE II: The numerical results for the scaling law. 



various a, N'j has been tabulated in the fifth column of 
the same table. In order to arrive at the scaling laws, we 
expand the left hand side of Eq. (|16p in powers of m? / h? 
and keep the leading terms. Thus 



-2/s 



V{A)T{A + 2) 



{a - 7^) 

2777 



-lB [Ja-i{B) - Ja+i{B)] - 2(1 - s - \)Ja{B) 
-iB[YA-i{B)-YA+i{B)]-2{l-s-\)YA{B)' ^ ' 

Carrying out a Taylor expansion near the critical cou- 
pling, we find the following scaling law 



m/A - K{Nf) (a-a,)i/2 



(19) 



where ha{Nf) = \/df/da\oi=a^- As anticipated, the 
scaling law for Nf comes out to be of the form : 



m/A = hNfia) {N^-Nff/^ , 



(20) 



witlr hNf{a) = ^df / dNf \nj=ni- These analytical re- 
sults are based upon the linearization of the original prob- 
lem and the identification of M [p^ — > 0) = m. Exact nu- 
merical analysis of the original non-linearized version of 




m/A= 0. 2888(1. -N,/1.0253)"'i2'' \i 
m/A= 0. 1813(1. -N,/l. 0253)° = 



FIG. 4: The scaling law for the dynamically generated mass 
as a function of Nf . The solid line is the fit to the numerical 
findings whereas the dashed line is the analytical result of the 
hnearized equation for the mass function. 

M(p^=0) /A vs a and Nf 




FIG. 5: The criticality plane for the dynamical mass in the 
phase space of a and Nf. The points shown are the numerical 
results obtained. 



Eq. (|10p confirms the qualitative nature of the aforemen- 
tioned analytical results. In Figs. ([T]) and ([2]), we depict 
the mass functions for different values of a and Nf. A 
study of the dependence of M(0) = m as a function of 
a and Nf, exemplified in Figs. ^ and ([3]), permits us 
to decipher the corresponding critical values and scaling 
laws. Tire numerical details differ slightly from the ana- 
lytical findings. We now have the scaling laws 

,71/A = a„(7V/)(a-a,)i/f^(^/) , (21) 
m/A = a^^.(a)(iV; -7V;)i/''"/(") . (22) 

Though none of ba{Nf) and 6jVy(a) is strictly 2 for a 



5 



broader fit of the scaling law (see Table HT]) . figures ^ 
and (|3|) show that the analytical results are not a bad rep- 
resentation of the exact results in the immediate vicin- 
ity of the critical coupling (compare the solid curves 
against the dashed ones). The infinite order phase tran- 
sition of the quenched QED softens out to a finite or- 
der transition in its unquenched version. It also has 
consequences for the mass anomalous dimensions. If 
B w C = — 1/4|, the large momentum behaviour 
of Eq. (fn|) implies 7,„ « 1 -|- s. On the other hand, if 
B >> C, 7m ~ 1 + s/2. The quenched limit trivially 
follows. The numerical analysis of the full equation also 
yields 7„i > 1. In order to obtain a finite electron mass in 
the limit of A — > oo, one requires charge renormalization 
(26j . Therefore, in this limit, we impose 



a(A) 



1 



(a^iV/))""^"^^ 
It implies the following /3-function 

P{a) = A da/dA = -ho,{Nf) [a - a^) 



(23) 



(24) 



Thus the /3-function has a stable zero at the point 
a = ttc, the result also obtained with the one loop 
approximation to the photon propagator, Q. Not 
only does a critical value of coupling separate chirally 
asymmetric and symmetric phases but so does also a 
critical number of flavours above which fermions cease 
to posses mass just like in QCD and QED3. Fig. ([5]) 
shows the criticality plane in the phase space of a and 
Nf. It interpolates and extrapolates the points obtained 
through the numerical analysis. We believe that our 
analysis can and should be extended to the study of 
QCD through its SDEs. This is for the future. 
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